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Abstract 


A  turbulent  mixing  layer,  presumed  to  be  caused  by  strong  shear  due  to  inertial  waves,  was 
observed  from  the  research  submarine  USS  Dolphin.  As  a  consequence  of  the  mixing,  a  density 
flux  was  set  up.  Although  inertial  waves  exhibit  motion  on  a  much  larger  scale  than  turbulence, 
the  horizontal  extent  of  the  waves  has  a  limit.  Therefore,  a  gravitational  imbalance  between 
the  turbulent  layer  and  its  outside  is  expected.  Data  indicated  that  the  gravitational  imbalance 
created  an  intrusion,  at  the  head  of  which  high  wavenumber  internal  waves  were  observed. 


1 


1.  Introduction 


Atmospheric  disturbances,  such  as  the  passage  of  a  storm,  in  the  open  ocean  often  cause  inertial 
waves  to  propagate  downward  (Price  1983;  Pollard  1970;  Pollard  and  Millard  1970;  D’Asaro 
1985,  1989).  Since  the  physical  scales  of  inertial  waves  are  at  least  several  kilometers  and  the 
total  amount  of  energy  is  very  high,  the  associated  shear  can  induce  a  persistent  turbulent  layer 
(Itsweire  et  al.  1989;  Hebert  and  Mourn  1994).  Usually  the  vertical  group  speed  of  near-inertial 
waves  is  an  order  of  several  meters  per  day,  thus  the  associated  disturbance  propagates  slowly 
in  the  water  column.  Although  inertial  waves  have  been  observed  at  all  levels  of  the  open  ocean 
water  column  (Sanford  1991),  turbulence  caused  by  near-inertial  waves  rapidly  dissipates  the 
kinetic  energy  of  the  wave  (Eriksen  1991;  Hebert  and  Mourn  1994). 

Regarding  turbulent  mixing  events  induced  by  near-inertial  waves,  we  need  to  bear  two 
facts  in  mind:  (1)  the  inertial  wave  has  a  limited  horizontal  scale  and  (2)  turbulence  elevates 
the  potential  energy  of  the  water  column.  Hence  we  can  expect  that  the  horizontally  linoited 
water  column  with  higher  potential  energy  than  the  surrounding  water  column  may  cause  a 
gravitational  collapse.  In  fact,  multisensor  observations  in  the  atmosphere  (Mahapatra  et  al. 
1991)  clearly  showed  an  undular  bore  generated  from  a  storm  mixing  event,  and  upstream 
disturbances  were  observed  in  front  of  the  bore.  A  similar  phenomenon  known  as  “the  morning 
glory”  has  also  been  reported  (e.g.  Smith  et  al.  1982).  If  the  same  scenario  is  applied  to  the 
ocean,  turbulent  mixing  is  not  the  end  of  the  energy  cascade.  The  vertical  density  flux  due  to 
turbulent  mixing  may  be  a  source  of  energy  for  a  gravitational  collapse  and  soliton-like  waves  can 
be  generated  from  the  upstream  disturbance,  as  has  been  observed  in  a  laboratory  (Maxworthy 
1980;  Amen  and  Maxworthy  1980).  This  process  can  generate  high  wavenumber  internal  waves 
away  from  the  solid  boundary  in  the  central  area  of  the  open  ocean. 

The  events  of  a  scenario  are  illustrated  in  Figure  1.  The  whole  process  is  as  follows: 

1.  an  atmospheric  disturbance  causes  an  inertial  wave; 

2.  turbulence  is  driven  by  the  inertial  wave  shear; 
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3.  vertical  density  fluxes  result  from  the  mixing; 


4.  a  gravitational  imbalance  in  the  horizontal  direction  induces  an  intrusion; 

5.  high  wavenumber  internal  waves  are  generated  at  the  head  of  the  intrusion. 

Making  use  of  data  published  by  Yamazaki  and  Osborn  (1993),  evidence  has  been  found 
to  justify  the  process  described  above.  The  next  section  briefly  summarizes  the  data,  with 
consistency  tests  presented  in  the  following  section. 

2,  Observed  Data 

Microstructure  data  were  observed  from  the  research  submarine  USS  Dolphin  which  has  a  length 
of  51  m  and  a  diameter  of  5.6  m.  Turbulence  was  measured  at  the  top  of  a  4.7-m  tripod 
placed  near  the  bow  of  the  submarine.  Two  airfoil  probes  were  used  to  measure  the  vertical 
component,  w,  and  the  horizontal  component,  u,  of  turbulent  velocity.  Two  high-response 
thermistors  measured  the  temperature  field  adjacent  to  the  airfoil  probes.  The  hydrography 
was  monitored  from  a  CTD  package  mounted  at  the  top  of  the  tripod.  Along  the  front  leg  of 
the  tripod,  six  thermistors  were  mounted  to  resolve  the  vertical  structure  of  stratification.  A 
1.2-MHz  ADCP,  with  a  1-m  resolution,  was  mounted  in  front  of  the  tripod  to  monitor  the  shear 
field.  Details  of  the  instrumentation  can  be  found  in  Yamazaki  and  Osborn  (1993). 

The  data  used  in  this  paper  come  from  part  of  a  7-h  dive  which  took  place  off  the  coast  of  San 
Diego,  California.  The  distance  to  the  shore  was  about  30  Km,  and  the  depth  generally  exceeded 
800m.  Two  ascending  and  two  descending  legs  of  the  tripod  located  between  the  surface  and 
about  50  m  in  depth  provided  the  vertical  structure  of  the  background  stratification.  Although 
the  dive  took  place  at  night,  the  effects  from  convective  cooling  was  limited  to  above  30  m  (Figure 
2).  During  0950  and  1030  UTC,  we  encountered  a  vigorous  shear  layer  with  the  shear  sufficiently 
strong  enough  to  induce  an  energetic  mixing  event  (Yamazaki  and  Osborn  1993).  The  turbulent 
shear  layer  had  a  thickness  of  about  5  m  and  was  located  in  a  seasonal  thermocline  (Figure 
2).  Because  the  salinity  contribution  to  the  density  was  minimal,  as  discussed  in  Yamazaki  and 


Osborn  (1993),  the  thermal  structure  obtained  from  a  thermistor  chain  provided  a  picture  of 
the  local  stratification. 

A  2400-s  data  segment  was  extracted  from  the  whole  dive.  Since  the  nominal  speed  of  the 
Dolphin  was  1.5  m  s~^,  the  displayed  segment  corresponds  to  3600  m.  The  thermal  structure 
was  monitored  from  the  thermistor  chain  (Figure  3a).  Turbulence  was  measured  along  the  solid 
line  at  the  top  of  the  temperature  contour.  Figure  3b  shows  the  high  shear  region  identified 
from  ADCP  data.  Again,  the  solid  line  in  this  figure  shows  the  location  where  turbulence  was 
observed.  Both  the  rate  of  dissipation  of  kinetic  energy,  e,  and  the  temperature  gradient  dissi¬ 
pation  rate,  x,  were  used  to  infer  the  mixing  layer  shown  in  this  figure.  The  local  hydrography 
conditions,  i.e.  the  local  gradient  Richardson  number,  Ri,  the  square  of  the  horizontal  shear,  Sh, 
the  local  buoyancy  frequency,  N,  salinity,  S,  density,  Cf,  and  temperature,  T,  were  computed 
from  the  CTD  data  set  and  10-s  average  ADCP  data.  See  details  of  the  data  processing  in 
Yamazaki  and  Osborn  (1993). 

3.  Consistency  Tests 

In  order  to  identify  the  source  of  the  horizontal  shear,  two  sections  of  the  shear  vector  field 
were  examined;  1)  between  the  elapsed  time  of  500  and  1000  s,  where  the  turbulence  probes 
were  in  the  mixing  layer  and  2)  between  the  elapsed  time  of  1200  and  1500  s,  where  the  shear 
magnitude  contour  is  the  most  contiguous.  The  shear  vector  measured  from  the  direction  of  the 
submarine  operation  showed  -120  at  the  top  of  the  layer  to  -80  toward  the  bottom  layer  for 
both  sections  (Figure  4).  This  shear  field  is  consistent  with  a  downward-propagating  inertial 
wave.  From  a  previous  experiment  in  the  same  area  (Yamazaki  and  Lueck  1987)  a  persistent 
shear  layer  was  fonnd  at  roughly  200  m  at  an  interface  between  the  California  current  and  the 
California  undercurrent.  No  shear  layer  due  to  any  other  type  of  flow  system  is  recognized 
above  100  m.  Therefore,  the  observed  shear  layer  is  due  to  a  near-inertial  wave  generated  from 
an  atmospheric  disturbance.  Similar  features  are  reported  in  Itsweire  et  al.  (1989)  and  Hebert 
and  Mourn  (1994).  If  this  inertial  wave  was  topographically  generated  from  a  tidal  forcing. 
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so  that  the  mixing  time  scale,  ^mix?  ^nay  be  estimated  from  the  following  discrete  form 


AT  __  {wT) 


where  AT  is  the  temperature  difference  between  the  top  and  the  bottom  of  the  mixing  layer,  and 
Az  is  half  of  the  mixing  layer  thickness.  The  observed  values  shown  in  Table  1  give  ^mix  =  106  h. 
Since  this  time  scale  is  solely  based  on  the  observed  vertical  heat  flux,  this  time  scale  may  be 
considered  to  be  the  upper  bound  for  the  mixing  time  scale. 

Because  the  edge  of  the  mixing  layer  was  observed,  we  cannot  assume  the  horizontal  homo¬ 
geneity  of  the  mixing  layer.  Therefore,  the  advection  terms  as  well  as  the  horizontal  flux  terms 
should  be  taken  into  consideration.  In  order  to  account  for  these  terms,  the  heat  flux  through 
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the  density  interface  at  the  bottom  of  the  mixed  layer  is  estimated.  The  reason  for  doing  so  is 
that  both  the  horizontal  terms  and  the  vertical  fluxes  within  the  mixed  layer  should  result  from 
the  heat  flux  through  the  interface.  An  attempt  to  calculate  the  heat  flux  is  made  using  the 
eddy  diffusivity  method  of  Osborn  and  Cox  (1972).  Since  the  vertical  flux  within  the  mixing 
layer  has  been  observed  (Yamazaki  and  Osborn  1993),  the  eddy  diffusivity  can  be  estimated 
from 


{wT') 

(dT„ldz)„ 

where  {dToldz)m  is  the  temperature  gradient  in  the  mixing  layer, 
interface  is  expressed  as 


(4) 


and  the  heat  flux  through  the 


(5) 


where  {dTofdz)i  is  the  temperature  gradient  at  the  bottom  of  the  mixing  layer  interface.  Here 
the  eddy  diffusivity  at  the  interface  is  assumed  to  be  identical  to  the  eddy  diffusivity  in  the 
mixing  layer.  The  mean  background  temperature  gradient  in  the  mixed  layer  is  0.03°C  m“^  so 
Kt  is  8  X  10“®  m^  s“^.  Making  use  of  this  coefficient,  the  heat  flux  at  the  density  interface  is 
expected  to  be  4.4  x  10~®  °C  m  s“^,  so  the  mixing  time  scale  becomes  3.9  h.  At  least  this  is 
the  minimum  time  scale  to  create  the  mixed  layer.  In  order  to  induce  gravitational  collapse, 
turbulent  mixing  must  persist  longer  than  this  time  scale.  Thus,  this  is  the  lower  bound  for  the 
mixing  time  scale.  Since  inertial  wave  energy,  in  general,  propagates  at  several  meters  per  day, 
the  true  mixing  time  scale  is  in  a  reasonable  range  specified  by  the  upper  and  the  lower  bound 
time  scale. 

Although  the  submarine  fluctuated  slightly  between  37  and  41  m,  both  inside  (0-1000  s) 
and  outside  (1600-2400  s)  of  the  mixing  layer  were  observed  at  roughly  the  same  depth.  The 
difference  in  at  between  the  inside  and  outside  of  the  mixing  layer  is  about  0.05,  so  the  water 
column  is  gravitationally  unstable  in  the  horizontal  direction.  Therefore,  it  is  not  unreasonable  to 
expect  a  gravitational  collapse  caused  by  the  turbulent  mixing  layer.  The  gravitational  collapse 
causes  an  intrusion  into  a  stratified  fluid.  The  intrusion  takes  place  from  left  to  right  in  Figure 
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3,  so  internal  waves  between  the  elapsed  time  of  1200  and  1500  s  are  associated  with  the  head 
of  the  intrusion. 

In  general,  the  intrusion  speed  cannot  exceed  a  critical  speed  of  gravity  current  Cq 


where  h  is  the  thickness  of  the  gravity  current,  Ap  is  the  difference  between  the  gravity  current 
and  the  ambient  water,  po  is  the  ambient  water  density,  and  g  is  the  gravitational  acceleration. 
So  the  intrusion  speed  Ci  should  be  less  than  Cq  —  0.05  m  s“^. 

Cl  <  Cg^  (7) 

The  observed  internal  waves  and  the  intrusion  are  consistent  with  the  laboratory  experiment 
of  Amen  and  Maxworthy  (1980)  who  found  a  series  of  solitary  waves  generated  from  the  head 
of  a  gravitationally  collapsed  mixed  layer.  The  largest  amplitude  internal  wave  had  an  apparent 
wave  length  of  roughly  La  =  70  m.  Since  we  cannot  infer  at  what  angle  the  submarine  crossed 
the  internal  waves,  the  actual  wave  length  should  be  less  than  the  apparent  one.  Making  use  of 
the  observed  stratification  at  this  point  of  the  water  column,  N  —  0.01  s“^,  the  maximum  phase 
speed  of  this  wave  is 

C'max  =  LaN/27r  =  0.11  m  s“\  (8) 

The  largest  observed  wave  height  is  roughly  3  m,  so  the  minimum  wave  length,  Xmin?  correspond¬ 
ing  to  the  steepest  wave  slope,  0.34/7r,  (Thorpe  1978)  is  roughly  30  m.  Hence  the  minimum 
phase  speed,  Cmin.  is  0.05  m  s”^.  The  actual  phase  speed,  C,  should  be  between  Cmin  ^nd  Cniax- 

C*min  <  (7  <  C^max*  (9) 

Therefore,  observed  internal  waves  are  propagating  faster  than  the  intrusion. 

Based  on  the  foregoing  consistency  test,  it  is  concluded  that  the  proposed  energy  cascade  is 
a  plausible  process  in  the  ocean.  The  implication  of  this  result  is  that  turbulence  mixing  can 
cause  a  secondary  process,  i.e.  the  intrusion  due  to  gravitational  collapse. 
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Table  1.  Observed  values  associated  with  the  mixing  layer. 


Parameter 

Value 

{w'T') 

1.63  X  10-®  °C  m  s-i 

AT 

0.25 

Ap 

0.05  kg  m-3 

h 

5  m 

{dToldz)m 

0.03  °C  m-i 

{dToldz)i 

0.55  °C  m-i 
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Figure  L  Illustration  depicting  the  energy  cascade  proposed  in  this  paper.  An  atmospheric 
disturbance  generates  near-inertial  waves  whose  group  speed  Cg  is  almost  vertical  and 
whose  frequency  is  roughly  local  /.  Shear  associated  with  the  inertial  wave  causes  turbulent 
mixing  and  induces  the  density  flux  {wp).  The  flux  reaches  a  gravitational  collapse  of 
the  mixed  layer;  as  a  result,  the  head  of  the  intrusion  causes  the  upstream  disturbance 
generating  high  wavenumber  internal  waves. 

Figure  2.  Temperature  profiles  obtained  from  the  second  descending  (thick  curve)  and  the 
second  ascending  (thin  curve)  sections  of  the  dive.  A  shear  layer  was  observed  between 
35  and  45  m  during  the  descending  section,  and  no  shear  layer  was  observed  during  the 
ascending  section. 

Figure  3,  Hydrographic  conditions  associated  with  the  turbulent  mixing  layer,  (a)  Thermal 
structures  identified  with  a  thermistor  chain.  The  temperature  contour  interval  is  0.05^C. 
Making  use  of  turbulent  dissipation  rate  data,  the  depicted  mixing  area  is  shown  with 
hatches,  (b)  Shear  contours.  Two  levels  of  shear  contour  are  shown:  Sjj  =  4  X  10“^  8““^ 
and  8  x  10“^  s“^.  The  shaded  area  depicts  the  high  shear  region.  Thick  solid  line  in 
(a)  and  (b)  is  the  location  of  the  turbulence  package,  (c)  The  dissipation  rate  of  thermal 
variance  x*  (d)  The  rate  of  dissipation  of  kinetic  energy,  c.  (e)  Local  gradient  Richardson 
number,  Ri.  (f)  Local  shear  squared  value,  (g)  Salinity,  S.  (h)  at.  (i)  Temperature, 
T.  (j)  Local  buoyancy  frequency,  N^. 

Figure  4.  Contours  of  the  direction  of  the  shear  vector  are  drawn  for  (a)  elapsed  time  between 
500  and  1000,  and  (b)  between  1200  and  1700.  The  angle  is  measured  in  the  horizontal 
plane,  relative  to  the  direction  of  submarine  travel,  as  shown  in  (c).  The  shear  vector 
rotates  clockwise  as  it  descends.  The  background  contours  in  (a)  and  (b)  show  the  mag¬ 
nitude  of  shear  squared  appearing  in  Figure  3b.  The  location  of  the  turbulence  package  is 
shown  with  dashed  line. 
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